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We calculate the single transverse-spin asymmetry for open charm production in pp collisions 
within the QCD coUinear factorization approach. We include contributions from both twist-three 
quark-gluon and tri-gluon correlation functions. We find that the quark-gluon correlation functions 
alone generate only a very small asymmetry for open charm production in the kinematic region of 
current interest at RHIC, so that the observation of any significant single-spin asymmetry would be 
a clear indication of the presence of tri-gluon correlations inside a polarized proton. We furthermore 
demonstrate that the tri-gluon contribution could be very different for the production of D and D 
mesons. These features make the single spin asymmetry in open charm production in polarized pp 
collisions at RHIC an excellent probe of tri-gluon correlation functions. 

PACS numbers: 12.38.Bx, 12. 39. St, 13.85. Ni, 14.40.Lb 



I. INTRODUCTION 



Large single transverse-spin asymmetries (SSAs) have been consistently observed in various experiments at different 
collision energies [U, 0, 0] . These experimental data have triggered much theoretical activity aiming at identifying the 
fundamental mechanism behind the measured asymmetries as well as the relevant components of hadron structure. 
Among the proposed theoretical mechanisms, the transverse momentum dependent (TMD) parton distribution ap- 
proach d, [^, 0, H, H, and the tvsfist-3 quark-gluon correlation approach [ll|, [l^, [H, have been the ones most 
widely discussed in the literature in recent years. The two approaches each have their kinematic domain of validity. 
For semi-inclusive deep inelastic scattering (SIDIS) and for the Drell-Yan process, they were shown to be consistent 
with each other in the kinematic regime where they both app ly Q. Both approaches have been applied extensively 
in phenomenological studies [H, \Mm, [H, [H, S El, ■ 

For single hadron inclusive production in hadronic collisions, the large transverse momentum Ph± of the produced 
hadron is the only hard momentum scale. For the SSAs in such single-scale hard processes, one expects the QCD 
collinear factorization approach to be valid [1^. It attributes the observed phenomenon of SSAs to three-parton 
correlations inside a polarized proton, commensurate with the twist-3 nature of the overall observable. At the leading 
power in the hard scale, such partonic correlations are represented by the twist-3 quark-gluon correlation function 

tni,[ii, 

and by a twist-3 tri-gluon correlation function [23, [2^ , 

Tg{x,x) = j ^^L^e"^"'^r_i_(P,,^|F+„(0) [e^-'^^'' F+ {y^)] F'^^ {y^)\P, st) , (2) 

where the proper gauge links have been suppressed. There is a quark-gluon correlation function, Tq^p, for each quark 
(anti-quark) flavor q (q), and there are two independent tri-gluon correlation functions, Tq\x,x) and Tq\x,x), 
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because of the fact that the color of the three gluon field strengths in Eq. ([2]) can be neutralized by contracting with 
either the antisymmetric if^^^ or the symmetric d"^^^ tensors with color indices, A, B, and C [IJ, UHl- 

Most phenomenological studies of SSAs in the twist-3 formalism performed so far have concentrated on the con- 
tribution by the quark-gluon correlation functions [H, [3, H^- The role of the tri-gluon correlation function in SSAs 
was first studied by Ji in the context of direct-photon production in hadronic collisions [23 |. Recently, two of us 
investigated the contribution of this correlation function to the SSA for open charm production in SIDIS [25l |. The 
twist-3 tri-gluon correlation functions have also been discussed in the context of spin asymmetries for di-jcts (or 
di-hadrons) in hadronic scattering, calculated in the TMD framework (26j . Here they appear after weighting with the 
transverse-momentum imbalance of the two jets. 

In this paper, we study the role of the tri-gluon correlation function in the SSA for open charm production in 
hadron-hadron collisions. We calculate both the quark-gluon and the tri-gluon contributions to the SSA, at leading 
order. We find that the quark-gluon correlation function, extracted from data on SSAs in pion production, leads to 
a very small SSA for open charm production in pp collisions at RHIC energy in almost all the accessible kinematic 
region. That is, any sizable SSAs observed for open charm production at RHIC will be a direct evidence of a non- 
vanishing tri-gluon correlation inside a polarized proton. Within a simple model for the Tq\x,x) and Tq\x,x), we 
find a fairly large SSA for open charm production at RHIC energies. The SSA tends to be maximal in the forward 
region of the polarized proton and becomes small or changes sign in the backward region. In addition, we find that 
our calculated SSAs for D and D mesons are the same if the tri-gluon correlation function Tq'^ vanishes. If Tq'^ is 
finite, the SSAs for D and D meson production could be very different and could in fact be used to disentangle the 
contributions from the two different tri-gluon correlation functions. With data on SSAs for open charm production 
becoming available from RHIC [l^l , we will be able to extract the tri-gluon correlation functions and to learn about 
the dynamics of quantum correlations of gluons inside a polarized proton. With additional data on the SSAs for other 
pp processes Q and perhaps in the future for SSAs for D and D meson production in SIDIS, a global analysis of QCD 
dynamics beyond leading twist may become possible in the not too distant future. 

We note that predictions for the SSA in D meson production at RHIC have also been obtained within the TMD 
approach [2^, assuming that a TMD factorization holds for the single-inclusive observable pp — > DX. The asymmetry 
was considered as a means to learn about the gluon Sivers function. In our opinion, however, the twist-3 approach 
provides the relevant and more appropriate framework here, as we described above. We emphasize again that our 
approach naturally involves two separate tri-gluon correlation functions Tq\x,x) and Tq\x,x), whereas only one 
Sivers function occurs in the study of [2^ (see also [1^). This is expected in general to lead to discernible differences 
also in phenomenological applications. To give one example, using a single gluon Sivers function, one will find the 
SSAs for D and D meson production to be equal, whereas in our approach they could potentially be very different, 
as we shall see below. We will return to these issues in our concluding remarks. 

We also note that the recent work considers the SSA in charm and anti-charm production at RHIC, through 
the twist-3 mechanism. In this paper, only the contribution by the quark-gluon correlation function was considered, 
and the fragmentation of the charm quark into a charmed meson was neglected. 

The rest of our paper is organized as follows. In Sec. |TT1 we present our calculation of the SSAs for open charm 
production in hadronic collisions in terms of the QCD coUinear factorization approach. In Sec. Illli we propose simple 
models for the tri-gluon correlation functions, and present our numerical estimates for the SSAs for open charm 
production in pp collisions at RHIC. Finally, we summarize our results in Sec. IIVI 



II. CALCULATION OF SINGLE TRANSVERSE-SPIN ASYMMETRY 



We consider inclusive single charm meson production in a scattering process between a polarized proton A of 
momentum P and transverse spin vector st and an unpolarized proton B of momentum P', 

A{P, st) + B{P') ^ h{PH) + X, (3) 

where h represents the observed open charm [D or D) meson with momentum Ph and mass irih- 
The spin-averaged and spin-dependent cross sections <7{Ph), Aa{Ph, st) are defined as 

a{Ph) = i [a{Ph,ST) + <7{Ph, -st)] , Aa{Ph,ST) = ^ HPh^ST) ~ <j{Ph, -st)] ■ (4) 
The single transverse-spin asymmetry An is defined as the ratio of Aa{Ph, st) and a{Ph), 

dAa{Pu,ST) /„ da[P,,) 
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for the single hadron differential cross sections. 

Tlie spin-averaged differential cross section for D meson production at large transverse momentum, Pii± > to/j, can 
be written in the following factorized form [sil, [s^l : 



da 



a,b 



dz 



dx' 



dx 



(6) 



where J2a b represents the sum over all parton flavors and 5 = {P + P')^ is the total collision energy squared. 4>a/Aix) 
and 4>b/BW) 3're the standard parton distribution functions, and Dc^h{z) is the fragmentation function for a charm 
quark c fragmenting into a D meson. We have neglected all dependence on the factorization and rcnormalization 
scales in ([5]). When Phi^ 3> mu-, a sizable amount of charm mesons could be produced from the fragmentation of a 
gluon or a light quark or anti-quark. Then, the sum ^ in Eq. ([6]) would need to be complemented by an additional 
sum over the fragmenting parton, running over all possible light parton flavors as well as the charm. However, wc do 
not expect a large fragmentation contribution from light partons to charm meson production at the current RHIC 
energies, and we will limit our calculations to the fragmentation of charm and anti-charm quarks in this paper. 

In Eq. ([6]) , iJ^^j, is a short-distance hard part for two partons of flavor a and h to produce a charm quark c. At the 
lowest order, it gets contributions from the light quark-antiquark annihilation and gluon-gluon fusion subprocesses, 
as sketched in Fig. [1] and is given by 



ttU 

qq — >c 



H, 



Nc 
1 

27Vc 



1 

tu 



Ncl_ 
Cf P 



Amis — 



Amis' 
tu 



(7) 



where s,t,u are defined at the partonic level as 

S^{xP + x'P'f, i ^ {xP - pcf - ml, u = {x' P' - pcf - ml, (8) 
with pc and mc the momentum and mass of the charm quark that fragments into the D meson, respectively. 






(a) 



(b) 



FIG. 1: Lowest order Feynman diagram for light quark-antiquark annihilation 
quarks. 



and for gluon-gluon fusion to a pair of heavy 



The transverse spin-dependent cross section also gets contributions from both the qq annihilation and gg fusion 
channels. Following the generalized factorization theorem [23j . the cross section in the qq annihilation channel can be 
expressed in terms of a general twist-3 quark-gluon correlation function [2^ , 

dAa{sT) cc J dzDc^h{z) J —4'q/B{x') J dxidx2Tg^F{xi,X2) ie'"''^"'^ lh-n^-^j-p-Hqq^cixi,X2,k±), (9) 
with the quark-gluon correlation function Tq,p{xi, X2) given by 

Tq,F{x,,X2) = J ^^L^e-^^+^r+'(-^--^)^+^M^,STK\(0)7+ [e^^'^"'' F,+ (y-)] V^,(y-)|P, s^) , (10) 

whose diagonal term, at xi = X2, is equal to the correlation function in Eq. ([T]). In order to generate a non- vanishing 
SSA, a strong interaction phase is necessary, which comes from the interference between a real part of the scattering 
amplitude and an imaginary part of the partonic scattering amplitude with an extra gluon, as shown in Fig. [2l 
Technically, the imaginary part arises when the virtual momentum integral of the extra gluon is evaluated by the 
residue of an unpinchcd pole from a propagator in the amplitude with an extra gluon. Such a propagator is indicated 
by the short bars in the diagrams in Fig. O The phase can arise from the attachment of the extra gluon to cither the 







FIG. 2: Feynman diagrams that give the twist-3 contribution to the spin-dependent cross section in the quark-antiquark 
annihilation channel: initial-state interaction (a), (b), and final-state interaction (c), (d). The short bar indicates the propagator 
that produces the unpinched pole. 

initial-state parton, or the final-state charm quark, which we will refer to as initial-state and final-state interactions, 
respectively. 

At lowest order, there are four diagrams contributing to the twist-3 polarized cross section in the quark-antiquark 
annihilation channel, as sketched in Fig.[2l where the blob is given by the diagram in Fig.[T]Ja). Using the techniques 
well established in literature (see, for example, [131), we have 



Ej 



dAa 



a 
5 



dz 



dx' 

— (f>q/Bix') 



dx 



^PhSTTin • 



Tq,F{x,x) - X — Tq,p{x,x) ) H, 



^ 5{s + t + u) 

qq^c{s, i, U) + Tq^F{x, x)Hqq^c{s, t, u) 



(11) 



where Ha 



can be written as 



Hqq^c Hqq^c 



Hqq^c ' 1 



(12) 



and likewise for T-lqq—,c, and where the corresponding hard parts are given by 



ttI 

qq — >c 



Hi 



'^qq—*c 



1 



2^c 



2mis 



2Nl 



2mi 



2N^ 
2Nl 



2mt 



(13) 



(14) 



Note that Tiqq^^ and TLqq^^ are proportional to the charm quark mass. As a check of our results, when 





q'q' 



the spin-dependent cross section in Eq. (jlip becomes identical to the one for pion production through the qq 
channel [22| (if one replaces the D meson fragmentation function by the pion fragmentation function). 

The spin-dependent cross section for D meson production can be calculated in the same way. The Feynman 
diagrams are the same as those for D meson production in Fig. [21 except that the extra gluon should be attached 
to the anti-charm c quark for the final-state interaction. The cross section for D meson production has the same 
factorized form as that in Eq. pT|) . with the fragmentation function Dc^n{z) replaced by Dg_,^(z), and the hard 
parts given by 



ttI 

qq — 



ttF 



^qq — >c 



1 



-Hi 



2N^ 



r + 2mis 



2ml 
s 

27711' 



(15) 



(16) 
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These short-distance hard parts are consistent with those presented in the calculation of the SSAs for heavy quark 
and anti-quark production in hadronic collisions [sot . We note that the hard parts for qq scattering are obtained from 



those for qq by hI'^^ = ^-^a^^e ^'^d H^^^^ = —H^^^^, and hkewise for the H^'^ . 



Similar to Eq. ([9]), the spin-dependent cross section for the gg fusion channel has the following factorized form: 



dAa{sT) oc 



1 

25 



dx' 



dzD,^h{z) / —4>g/B{x') / dxidx2TG{xi,X2) ie"^^"" lim 



dk 



H, 



(xi,X2,fc_L), (17) 



where Tg{xi, X2) is defined as 

P+dy^dy^ 



Tg{xi,X2) = 



27r 



(18) 



with daf} = —gap + nafip + npUa- Tg(xi,X2) is related to the tri-gluon correlation function through Tg{x,x) 
xTg{x, x). 







FIG. 3: Feynman diagrams that give the twist-3 contribution to the spin-dependent cross section in the gluon-gluon fusion 
channel: initial-state interaction (a), (b), and final-state interaction (c), (d). The short bar indicates the propagator that 
produces the pole. 



To calculate the partonic hard part, Hgg—,c, in Eq. (|17p . we need to consider Feynman diagrams with either initial- 
state or final-state interactions, as sketched in Fig. [31 where the blob is given by the sum of the three diagrams in 
Fig. [TJ^fe). Hence, each diagram in Fig. [3] corresponds to nine diagrams. Instead of four diagrams in Fig. [2] for the 
quark-antiquark annihilation subprocess. we have a total of 36 diagrams for gluon-gluon fusion. By evaluating these 
diagrams, we obtain the contribution to the spin-dependent cross section. 



E, 



dAa 



'd^Ph 



S 



E 

i=f.d 



dz 



dx' 



Dc^h{z) / —(t>g/B{x') / — V4ia 



dx 



.PhSTnn 



S [s + i+u) 



(t^^ {x,x)- x^^ (x, x)^ H^'^-.ciS, i, u) + {x, a:)H«-.c (S, I u) 



(19) 



where the sum, X^i^/ d' over the two correlation hmctions Tq\x, x) and Tq\x, x). The partonic hard part Hgg-,c 



can be written as 



and likewise for Tigg^c, and we find 



aa—*c ^^gg-^c I ^-^gg^c 1 



(20) 



yy^" 8Cf iuS^ 
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1 u-t 
8Cf tus 



r +U'' + Amis 
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tu 



tu 
1 



1 



and 



.4 (7V2 - 1) & ANc [Nl - 1) iiL 
1 ml (P + u'^) (iu - 2mls) 



P + iP + Amis - 



Amis 



4-2 



tu 



(21) 
(22) 



2Cf 
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- ' 2Cf 



t^i 
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(23) 
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The gluon-gluon subprocess of course also contributes to the cross section for D meson production. The corre- 
sponding partonic hard parts for producing an anti-charm quark are given by 

rr(/) = /7-(/) flid) ^ _ljid) 

gg-^c ^^gg^cJ ^^gg-^c ^^gg^cJ 

_ _ -zjif) jj{d) _ ^ _n,(d) /24) 

^gg^c ' ^gg~tc ; ' '■gg^c ' '■gg—tc ; V^^y 

where the sign difference of the partonic hard parts for the Tq'^ contribution will be responsible for the difference of 
the SSAs for D and D meson production that will be discussed in the next section. We note that this sign difference 
can also be observed in the expressions for the "gluonic pole matrix elements" given in [2^. 

We point out that the compact dependence of the spin-dependent cross section on the combinations T{x, x) — 
xT'{x, x) of the twist-3 correlation functions found in Ref. [2^ for the "masslcss" case of pion production in hadronic 
collisions, is violated for the production of D (or D) mesons by the additional non-derivative terms in Eqs. (jlip and 
pop . The violation is caused by the heavy quark mass since the additional terms vanish when ttIc ^ 0. In fact, we 
observe that the hard parts we have derived satisfy the following relation: 

separately for any of the various contributions considered above (and likewise for c production) . This connection is 
likely a consequence of the "master formula" for twist-3 soft-gluon-pole contributions derived in [T3 | . 

We also note that "soft-fermion pole" contributions [l^l , for which the pole in the hard-scattering function is taken 
in such a way that the initial quark, rather than the initial gluon, becomes soft, are absent for the qq process at 
the leading order. This is because qq annihilation proceeds through an s-channcl diagram, whereas soft-fermion 
poles would only appear in t-channel diagrams. If they were present, such contributions would involve the function 
Tq ^iO, x). For the tri-gluon correlation contribution, terms proportional to Tq \o, x) and Trf\o, x) are automatically 
included in our calculations, thanks to the symmetry properties of the correlation functions |24{ . 

Combining the factorized cross sections in Eqs. (jlip and (jl9p with the corresponding partonic hard parts, we have 



dAa dAa 



qq^cc " 



(26) 

gg^cc 



for the leading-order contribution to the transverse-spin-dependent cross section for D (or D) meson production in 
hadronic collisions. The corresponding single transverse-spin asymmetry is obtained by substituting Eqs. ^ and 
into Eq. jS]). 



III. PHENOMENOLOGY 



We now use the expressions we have derived in the previous section to obtain phenomenological results suitable 
for RHIC. We will propose simple models for the tri-gluon correlation functions Tq\x,x) and Tq\x,x), and then 
estimate the size of the SSA for D and D meson production in p^p collisions at RHIC at ^/s = 200 GeV. 

In our studies, we adopt the CTEQ6L parton distribution functions [s^ for the unpolarizcd proton. We use the 
charm-to- 1? meson fragmentation function of Ref. [s^ . In order to treat the kinematic charm mass effects in the 
fragmentation process, we adopt one of the choices introduced in Ref. [sst . which corresponds to setting P/^x = ~Pc± 
and yu = Tjc = y, where Ph_L iPc±) and yo (yc) are the transverse momentum and rapidity of the D meson (charm 
quark), respectively. With this choice, we then have for s,t,u and the Feynman variable xp- 

s = x'xS, t = ~xrar_\\fSe^'^ ^ u = ~x'mc±y/Se^ , xp = — j=- (e^ — e~^) , (27) 

where m^x — \/ + -Pc± f^h^ = \l nif^ + P^^ . We further assume that D^_^jj{z) ~ Dc^£,{z) for the D meson 
fragmentation functions. We choose the factorization scale to be equal to the renormalization scale throughout, and 
set jjL = \/m1+ P^^ with TOc = 1-3 GeV. 

In order to calculate numerical values for the SSAs for the production of D or D mesons, we need the unknown, but 
universal, tri-gluon correlation hmctions Tq''^\x,x). Because of their nonperturbative nature, both the quark-gluon 
correlation functions Tq^p{x,x) and the tri-gluon correlation functions Tq''^\x,x) should in principle be extracted 
from experimental data on SSAs. For the Tg^p{x, x), this was done in Ref. [1^ by fits to the experimental data [l|, Q 
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for the SSA in pp nX. We adopt in the following the set referred to as "Fit IF in [13] ■ We note that the tri-gluon 
correlations were not taken into account in p^ . so that no direct information on these is available so far. For the 
purpose of estimating the SSAs and motivating future experimental measurements, we therefore follow the arguments 
given in Ref. [1^ and model the T^''^\x, x) as 

tP{x,x) - XfG{x), T^g\x,x) = \aG{x), (28) 

with G{x) the ordinary unpolarizcd gluon distribution function. In order to cover a range of possibilities for the 
nonperturbative correlation functions, we introduce three sets of values for the parameters A/ and A^: (1) A/ = 

Ad = 0.07 GeV, (2) A/ = A^ = 0, and (3) A/ = -Ad = 0.07 GeV, corresponding to the assumptions: T^^ = T^f^, 

Tq^ = Tq^ = 0, and Tq'' = ~Tq \ respectively. In principle, the signs and the values of A/ and Ad, as well as 
the functional form of the correlation functions should be fixed by future data. Comparing SSAs for physical cross 
sections involving the same non-pcrturbative twist-three correlation functions but different partonic hard subprocesses 
will provide stringent tests of QCD dynamics and the twist-3 factorization we use. 




y y 

FIG. 4: The SSA as a function of rapidity y for meson (left) and _D° meson production (right) at ^/s = 200 GeV and 
Phx = 2 GeV. The curves are: solid (A/ = \d = 0.07 GeV), dashed (A/ = Ad = 0), dotted (A/ = -Ad = 0.07 GeV). 

In Figs. [4] and [5] we plot the SSAs, Ajv, for the production of D and D mesons as functions of rapidity y and 
Feynman-a;i?, respectively. We count positive rapidity in the forward direction of the polarized proton. The solid, 
dashed, and dotted curves correspond to the three sets of parameters: A/ = Ad = 0.07 GeV, A/ = Ad = 0, and 
\f — —Ad — 0.07 GeV, respectively. From the dashed curves in Figs. [4] and [5l it is clear that the quark-gluon 
correlation function Tq^p alone generates a very small single transverse-spin asymmetry at RHIC energy. This is 
because of the dominance of the gg fusion contribution over the qq one in the spin-averaged cross section in the 
denominator of Aiq. In other words, any significant size of the SSA in open charm production signals the discovery 
of tri-gluon correlations inside a polarized proton. 




FIG. 5: Same as Fig. U but as a function of Feynman-XF. 
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The difference between the solid and dotted curves in Figs. |4] and [5] indicates that the two tri-gluon correlation 



functions, Tq and Tq' , may both play very important, but different, roles for the SSA in D and D meson production. 

In the case of D mesons, a large (see the solid curve) is obtained when Xf = Xd = 0.07 GeV, i.e., when Tq 

and have the same sign. However, when their signs are opposite, their contributions to the SSA tend to cancel, 
leading to a much smaller SSA (dotted curve). On the contrary, for D meson production, the largest A]\j is found 

when Tq'' and Tq"* have opposite signs. This is due to the fact that, as shown in Eq. ([24|) . the partonic hard parts 

associated with Tq*' change sign when going from charm to anti-charm production, while the hard parts for Tq^' 

remain the same. As a result, the SSA for D mesons is much smaller if Tq^ and Tq^ have the same sign. In addition, 
as seen from Figs. |4]and[5l scans of the SSA from the forward to the backward region may provide good sources of 
information on the x-depcndence or the functional form of the tri-gluon correlation functions, in particular if a sign 
change occurs. It is also striking to see that the asymmetry for either D or D mesons may become very large at 
forward rapidities at RHIC (but not for both simultaneously). 



(/) 




FIG. 6: The SSA as a function of Ph± for D" (left) and D° mesons (right) at mid-rapidity, y = 0, and 
curves are: solid (A/ = A<i = 0.07 GeV), dashed (A/ = = 0), dotted (A/ = - Ad = 0.07 GeV). 



200 GeV. The 



In Figs. [5] and [71 we show An for Z?" and l)" meson production as a function of Ph±, at mid-rapidity {y = 0) 
and forward-rapidity (y = 1.8), respectively. The absolute values of the SSAs decrease as a function of Ph±, which 
is a natural behavior of the twist-3 effect in QCD coUinear factorization. As before, while the contribution by the 
quark-gluon correlation functions is very small, the two tri-gluon correlation functions can make sizable, and very 
different, contributions to the SSA, thanks to the difference in the partonic hard parts in Eq. 



:0.08 



0.08 




FIG. 7: Same as Fig. |6l but at forward rapidity, y = 1.8. 



We finally note that first experimental data on the SSAs for open-charm production are now emerging from RHIC 
[2^ . The uncertainties are currently too large to allow distinction of the various models we have proposed, but the 
measurements are certainly very encouraging. 



9 



IV. SUMMARY AND CONCLUSIONS 



We have studied the single transverse-spin asymmetries for D and D meson production in hadronic coUisions. 
Within the QCD colhnear factorization approach, we have calculated both the derivative and the non-derivative 
contributions to the SSAs. We have included the contributions by the twist-3 quark-gluon correlation functions, 
shown in Eq. (fTTj) . as well as by the twist-3 tri-gluon correlation functions, given in Eq. ([19]). We have found that 
the quark-gluon correlation functions alone generate a negligible SSA for open charm production in almost all the 
kinematic region of interest at RHIC. We thus conclude that any sizable SSAs observed in D ot D production will be 
a discovery of tri-gluon correlations inside a polarized proton. 

We have proposed several simple models for the two unknown tri-gluon correlation functions and used these to make 
predictions for the SSAs in D meson and D meson production in polarized proton-proton collisions at RHIC. We 
have found that within our models, the asymmetries could be sizable and provide an important source of information 
on the two unknown tri-gluon correlation functions. We also found that the two correlation functions play a very 
different role in generating SSAs for D and D mesons, due to the sign difference between the partonic hard parts for 
charm and anti-charm production (see Eq. ([M)) ). Comparison of the SSAs for D and D mesons could thus provide an 
excellent tool for separating the two tri-gluon correlation functions. With data for the SSAs in open charm production 



arriving from RHIC [27[, we will be able to learn for the first time about the dynamics of multi-gluon correlations 
inside a polarized proton. 

We close this paper with a theoretical observation. As we mentioned in the Introduction, two approaches are 
commonly used for the study of SSAs in high energy collisions: the colhnear factorization we have used in this paper, 
and the transverse momentum dependent (TMD) factorization approach. The two arc known to be closely connected 
and complementary to each other 0]. As was shown in Rcf. Q, the twist-three quark-gluon correlation function, 
Tq^p{x, x), is related to a moment in transverse momentum of the corresponding quark Sivers fimction, qxix, k±^): 

TqAx.x) = / <m^klqT{x,k^), (29) 
Mp J 

up to ultraviolet renormalization. Here Mp is a hadronic mass scale. Applying the same method used to derive 
the above relation to the TMD gluon distribution, we would obtain an identical relation between the tri-gluon 
correlation function Tq\x,x) and the moment of the gluon Sivers function. The occurrence of the antisymmetric 
color contraction of the three gluonic field strengths in the tri-gluon correlation function is a natural consequence 
of the expansion of the gauge link in the adjoint representation when taking the transverse momentum moment of 
the gluon Sivers function. It is interesting to note that the other tri-gluon correlation function, Tq\x,x), does not 
appear to have a simple relation connecting it to the operator definition of the gluonic Sivers function. That said, in 
the TMD approach one can define tri-gluon correlators with both the / and d color structures [l^ when considering 
general QCD hard-scattering processes, which correspond to our Tq"^ and T^"^ after transverse-momentum weighting. 

In any case, both T^^\x, x) and T^^\x, x) are allowed and in general present in the colhnear factorization approach. 
It is therefore very important to measure the difference of SSAs in D and D meson production, in order to identify 
the role ofT^^\x,x). 
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